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CN Abstract 

^-1' Detection of electromagnetic (EM) counterparts of pre-coalescence binaries has very important implications for our understand- 
ing of the evolution of these systems as well as the associated accretion physics. In addition, a combination of EM and gravitational 
wave signatures observed from coalescing supermassive black hole binaries (SBHBs) would provide independent measurements of 
redshift and luminosity distance, thus allowing for high precision cosmological measurements. However, a statistically significant 
I— I sample of these objects is yet to be attained and finding them observationally has proven to be a difficult task. Here we discuss ex- 
isting observational evidence and how further advancements in the theoretical understanding of observational signatures of SBHBs 
before and after the coalescence can help in future searches. 
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■1. Introduction 

The most direct evidence for the formation of supermas- 
sive black hole binaries (SBHBs) currently available through 
electromagnetic (EM) observations is the existence of black 
hole pairs spatially resolved on the sky. An example of such 
a system is NGC 6240, an ultra-luminous infrared galaxy in 
which the Chandra X-ray observatory revealed a merging pair 
of X-ray emitting act ive nuclei with a separation of ~1 kpc 
llKomossa et al. , 2003 ). The large observed separation implies 
that the two supermassive black holes (SBHs) are a pair of unre- 
lated massive objects orbiting in the potential of the host galaxy 
that are yet to reach the binary stage, in which the two SBHs 
will be bound by their mutual gravitational force. Practical 
.obstacles in the direct identification of close binaries via EM 
observations arise from the need for a very high spatial resolu- 
tion and accuracy in position measurements. Because of the 
small orbital separation of sub-parsec binaries and the small 
separation of recoiling SBHs from the centers of their host 
galaxies not too long after the coalescenc^l the best prospect 
for direct imaging of these objects is offered by radio interfer- 
ometers, provided that SBHBs and recoiling SBHs are radio 
sources. Since this method relies on pointed observations of 
candidate objects, it requires a priori knowledge about their 
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location on the sky. This technique yielded one binary can- 
didate, 0402-1-379, an object with two compact, unresolved 
radio-cores at a projected separation of o nly ~ 7 pc, observed 



with the Very Lo ng Baseline Array (VLBA; Maness et al.L 12004 



Rodriguez et all 12006) . Followup multi-wavelength observa- 



tions of this and similar radio candidates in the future are essen- 
tial for studying the rest of associated observational signatures. 

The second most compelling piece of observational evidence 
for a SBHB is a periodic signal associated with Keplerian mo- 
tion, such as the periodic or quasi-periodic variability of the 
light emitted by an object. A binary candidate selected accord- 
ing to this criterion is the blazar OJ 287, which exhibits outburst 
activity in its optical l ight curve with a period close to 12 years 



(iValtonen et al. 



20081) . Distinguishing periodic variability due 



to a binary from the common variability of active galactic nuclei 
(AGNs) can be challenging because it relies on long-term mon- 
itoring of the candidates and is, in many cases, less pronounced 
than in OJ 287. 

The third technique used to select SBHB candidates relies 
on a detection of the Doppler-shift in the spectrum of an object 
that arises from the orbital motion of a binary. In this case the 
search utilizes the emission-lines associated with multiple ve- 
locity systems in the spectrum of a candidate object. Because 
this effect is also expected to arise in case of a recoiling SBH re- 
ceding from its host galaxy, the same approach is used to "flag" 
candidates of that type. The Doppler shift signature is how- 
ever not unique to these two physical scenarios, and comple- 
mentary observations are needed in order to determine the true 
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Figure 1: Illustration of the binary model for J0927. The innermost region of 
a circumbinary disk is shown after the binary has cleared a low density "hole" 
in the disk (top view, not drawn to scale). In this model r-NELs are associated 
with the circumbinary disk, b-BELs with the disk surrounding the less massive 
secondary SBH, and b-NELs with the accretion streams flowing from the in- 
ner edge of the ci rcumbinary disk toward the secondary. Figure adapted from 
iBogdanovic et alj |2009). 



nature of observed candidates (as discussed in 91.11 1. The ad- 
vantage of the method is its simplicity, as emission spectra that 
exhibit Doppler shift signatures are relatively straightforward 
to select from large archival data sets, such as the Sloan Digital 
Sky Survey (SDSS). So far a few SBHB candidates have been 
selected out of ~ 21, OOP SDSS quas a rs: J092712.65+2 9 4344 
(hereafter J0927; Bogdanovicetal., 2009; Dotti et al., 2009h. 
J153636.22+044127.0 (J1536: iBoros on & Lauer. 2009.) . and 
J105041. 35+34563 1.3 (J1050; ShieldseLaLll2Q09b|), where the 
first and the last object have also been flagged as recoi ling SBH 



candidates dKomossa et al.Ll2008t [Shields et al.L 12009b) . In the 



following section we describe in more detail results and un- 
certainties associated with this approach, using J0927 as an il- 
lustrative example, and summarize the status of the other two 
candidates. 

1.1. SBHB candidates: J0927, J 1536, and J 1050 

The optical spectrum of J0927 features two sets of emission 
lines offset by 2650 km s"' with respect to each other. The 
"redward" system consists of only narrow emission lines (r- 
NELs), while the "blueward" emission-line system comprises 
broad Balmer lines (b-BEL) and narrow, high-ionization for- 
bidden lines (b-NELs). In the context of the binary model, the 
observed velocity shift represents the projected orbital veloc- 
ity of the less massive member of a bound black hole pair. In 
this model, the emission lines of the blueward system (b-NELs 
and b-BELs) originate in gas associated with the less massive, 
secondary black hole, while the r-NELs originate in the inter- 
stellar matter (ISM) of the host galaxy. An illustration of the 
proposed geometry for the binary model is shown in Figure 1 . 
For an assumed inclination and orbital phase angles of 45° and 
a mass ratio of q - 0.1, the inferred orbital period and mass of 
the binary are ~ 190 yr and ~ 10^ Mq, respectively. The ac- 
creting secondary is the main source of ionizing radiation while 




Figure 2: Illus tration of the model of a recoiling SBH proposed for J0927 by 
Komossa et al. 12008) (oblique view, not drawn to scale). In this model r-NELs 
are associated with the ISM of the host galaxy, b-BELs with the disk recoiling 
with the SBH, and b-NELs with the expanding accretion disk, swept up ISM, 
or associated outflows. 



the primary is either quiescent or much fainter. 

According to the recoiling SBH interpretation for J0927 
(n ote that this was the fir st explanation proposed for this object 
by iKomossa et al.l 120081) . the b-BELs originate in the broad- 



line region retained by the SBH, the b-NELs were attributed to 
gas that is marginally bound to the SBH, swept up or outflow- 
ing, and the r-NELs are associated with the ISM of the host 
galaxy (the geometry of the recoiling SBH model is outlined in 
Figure 2). In this model the recoiling SBH of mass ~ 6x10^ Mq 
appears as an AGN during the period in which it accretes from 
the disk that is carried along with the hole. The accreting SBH 
provides a source of ionizing radiation for the b-BELs and b- 
NELs, as well as for the r-NELs, albeit from a larger distance. 

Another class of proposed explanations for J0927 includes 
interact ion of two galaxies in the potential well of a massive 
cluster (iHeckman e t al., 2009) and a super position of otherwis e 
unrelated AGN in a cluster environment dShields et al.Ll2009ah . 
Both of these scenarios can, in principle, explain the multiple 
velocity-emission-line systems seen in the spectrum of J0927 
and offer a good alternative to the SBH B and recoiling SBH 
models. In a subsequent paper, however. iDecarli et al. I (l2009ah 
reported that J0927 is unlikely to be a member of a rich cluster. 
At this point every model proposed for J0927 received some 
amount of criticism (see the above papers for details) and thus, 
the nature of J0927 is not conclusively determined and requires 
further investigation. 

Soon after its discovery, the interpretation for the candi- 
date J1536 was called into question because o f the lack of the 



variabihty expected f rom Keplerian motion dChornock et al 



2009al: iGaskel l2009h and the si milarity of its line profiles to 
those of "double-peaked emitters" (IChornock et all l2009b). At 
this point, the data available for this object do not conclusively 
rule out the binary scenario and further observational tests 
have been suggested dWrobel & Laon l2009t iLauer & Boroson , 



20091; iDecarliet al.ll2009bl ). 
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Most recently, 'Shields et al. ' ('2009b') reported the discovery 
of another unusual SDSS quasar, J1050, for which an SBHB is 
a viable explanation, along with a recoiUng SBH, and a double- 
peaked emitter. As with the previous cases, further observations 
are needed in order to discriminate among proposed hypothe- 
ses. 

J0927, J1536, and J1050 do not make up a uniform class 
of objects because their emission-line spectra differ from each 
other A common property shared by all three objects is that 
their emission line systems exhibit shifts in the range 2000 - 
4000 kms If such velocity shifts are interpreted in the con- 
text of a binary model, and if the emission-line systems are as- 
sociated with one or both of the black holes, the expectation is 
that the orbital motion of the binary should give rise to a mea- 
surable velocity change of the emission lines on a time-scale of 
several years. For all three objects this basic test was carried 
out by taking spectra of the object at a few different epochs and 
measuring the change in the position of the emission line peaks. 
In all three cases the measured rate of velocity shift is very low 
(dv/dt ~ few X lOkms"' yr"') and consistent with zero within 
the error bars. Taken at the face value this result seems to elim- 
inate the possibility that any of the three candidate objects is a 
binary. However, it may also be that if the selection process is 
biased towards binaries with large velocity shifts, they may be 
easiest to spot precisely when they are close to quadrature and 
dv/dt is small. Hence, it may take longer to detect a change in 
dv than expected if time scale is estimated based on a random 
phase prior Moreover, the uncertainties in emission geome- 
try and radiative processes of binary accretion flows currently 
preclude elimination of the binary model based on the lack of 
velocity shifts measured from the line peaks (see discussion in 

In the following section we briefly describe theoretical ideas 
for formation and evolution of SBHBs, including the conse- 
quences of their coalescence. 

2. Astrophysics of SBHBs before and after coalescence 

Galactic mergers are expected to be the major route 
for formation of binary black holes, in agreement with 
predictions of hierarchical model s of str ucture formation 
faaehnelt & Kauffmanni I2OO2I: IVolonteri et all l2003h and 
the observation that the majority of galaxies harbor massive 
black holes in their c e nters (e. g., [Ko rmendv & Richstone, 



for three-body interactions (IBerczik et al. ■ l2006t ISesana et al 
200 7a; Perets et a l.. 20071) and the availabiUty of gas which 
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Dotti et al., 20061 l2007l Maver et al., 20071 ICoM et al.. 


2007; 


Cuadraetal.1 200' 


i HavasakiL 2009h. 


In the third and final 



stage (at separations < 10"^ pc), the binary is led to coalesce by 
the emission of gravitational waves (GWs) and forms a single, 
remnant SBH with mass slightly lower (typically by several 
percent) than that of the two parent SBHs combined. Due to 
asymmetries in the orientation of the black hole spin axes with 
respect to the binary orbital axis, and the masses of the two 
SBHs with respect to each other, the emission of GWs is, in 
general, not isotropic. The GWs carry net linear momentum 
in some direction, causing the center of mass of the binary to 
recoil in the opposite direction. 

If the candidates described in ^1.11 turn out to be SBHBs, 
these would most likely correspond to binaries in the second 
evolutionary stage. Several physical processes may contribute 
to the presence of a, hopefully unique, spectral signature of 
such binaries: the presence of one or two (unresolved) accre- 
tion disks around individual black holes, albeit in general of 
different luminosity; the tidal distortion of the inner parts of 
the circumbinary disk by the presence of the binary (similar to 
features seen in proto-planetary disks). 

The "mass loss" and gravitational rocket effect at the end of 
the third evolutionary stage of the binary can have profound ef- 
fects on the surrounding gas and potentially give rise to unique 



observatio nal signatures as was recently suggested by a number 
of authors ("Loeb, '2007'; 'Shields & Bonning, 2008'; 'Lippai et al. , 
2008; Schnittman & Krolik, 2008; Kocsis& Loeb, 2008; 
Mo hayaee et all l2008t iDevecchi et all l2009t lO'Neillet al . 
I2009r Namely, the results of numerical relativity calculations 
show that recoil velocity of a remnant SBHs can range from 
about 200 km s ' for S BHs with low sp i ns or sp in axes aligned 



Baker et al 



about /UU km s tor a cHs with low sp i ns or sp m 
with t he orb ital axis ("Herrmann et alj, 2007allbt 
I2OO6I I2OO7I: [Gonzalez et al., 2007) to a remarkable -4000 
kms ' in case of maximally spinning SB Hs with spins op- 



19951: iRichstone et al.l Il998l: (Peterson & Wandell bOOOt 



positely directed and in the orbital plane (ICampanelli et al 
2007a b). If following a galactic merger the two SBHs find 
themselves in a gas rich environment, gas accretion torques 
will act to align their spin axes with the orbital axis, thus re- 
ducing the maximum kick to a value well below the galactic 
escape speed (the escape s peed from most galaxies is <2000 



Ferrarese & Fordi " 2005 ). Following a galactic merger, the kms Merritt et al. , 2004 ). This effect is expected to increase 



two SBHs are initially carried within their host bulges 
and eventually find th emselves at ~kpc separations (e.g.. 



the chance of reten t ion of recoiling SBHs by their host galaxies 



Kazantzidis et al.L 120051) from each other In this first evolu- 
tionary stage, the separation of a pair of the SBHs decreases 
through the process of dynamical friction that involves 
two-body gravitational interaction of individual SBHs with 
stars and interactions with gas in the nuclear region. In the 
second stage (at separations <10 pc), the binary becomes 
gravitationally bound when the binary mass exceeds that of the 
gas and stars enclosed within its orbit. The evolution of the 
binary in this phase is determined by the availability of stars 



(IBogdanovic et al.L l2007h . in agreement with the observation 
that a knost all gala xies with bulges appear to have a central 
SBH (IFerrarese & F ord. 2005). However, there may be some 
fraction of remnant SBHs formed in the aftermath of gas poor 
mergers, where gas torques are insufficient to prevent the high 
velocity recoil, launching the SBH out of the potential well of 
a host bulge or a host galaxy. One unique spectral signature 
of such systems could result from the distortion of a nuclear 
accretion disk and due to the velocity shift of the portion of 
the disk that is gravitationally bound to the SBH. On the basis 
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Figure 3: Snapshot from a simulation sliowing a SBHB and gas projected onto 
tlie plane of the binary orbit in the co-rotating model at 9.4 yr after the beginning 
of the simulation. The rotation of the binary and the disk is counter-clockwise. 
The temperatures of gas particles are marked with color: red T < 10"* K; yellow 
10-*K < T < 10* K; green lO^K < T < lO** K; and blue 10'*K <T < lO'" K. 
The higher temperature pailicles are plotted over the lower temperature ones, 
with the result that some information is hidden. Figure 4 ill ustrates this effect 
and can also be used as a color bar. Figure adapted from iBogdanovic et alj 
120081) . 



of this signature. iKomossa et alJ ( 120081) reported a discovery of 
the first re coiling SBH candidate, J0927, followed by J1050, 
reported by lShields etai.i C2009b) . (see ^TTTT ). 

Therefore, the emission-hne Doppler shift has so far been a 
key signature, sought after in screenings for both the SBHB and 
recoiling SBH candidates. In order to utilize it to its full poten- 
tial, and in anticipation of even reacher data sets in the future, 
expanded towards higher redshifts and lower-luminosity AGN, 
it is necessary to gain understanding of the remainder of associ- 
ated emission signatures in other wavelength bands. An oppor- 
tunity to study observational signatures in parallel with obser- 
vations is offered by hydrodynamics simulations of SBHs inter- 
acting with gas. In the following section we describe some of 
the advancements and challenges associated with this approach 
in general and then present selected results from our work. 

3. Simulations of SBHs and gas 

Over the past several years simulations of merging galaxies 
have significantly contributed to our understanding of the for- 
mation and evolution of black hole pairs (as discussed at the 
beginning of ®. However, simulations that would span with- 
out gaps the entire dynamical range from a galactic merger to 
binary coalescence are still computationally prohibitively ex- 
pensive. In order to reduce this problem into one that is nu- 
merically tractable it is currently necessary to make simplifica- 
tions in the initial conditions and the treatment of the thermo- 
dynamics of the gas. These choices include black hole pairs 
initialized in steady-state galactic and nuclear accretion disks, 
while in reality they are most likely to be perturbed and pos- 
sibly heavily irradiated by one or two AGNs. The treatment 
of gas thermodynamics in such disks may involve a parametric 



Figure 4: Left: Temperatures of gas particles as a function of radius at 9.4 
years in the co-rotating model simulation. This figure shows the extent of tem- 
perature stratification in the disk. Different colors that mark each temperature 
band are equivalent to those in the previous figure. This figure coiresponds 
to the morphology of the disk shown in Figure 3. Right: Histogram showirm 
the temperature distribution of particles. Figure adapted from Bog danovic et aU 
OOOS). 



prescription in the form of a fixed equation of state or a cool- 
ing law parametrized in terms of the dynamical time scale of a 
system (y6-model). These simplifying assumptions result in un- 
certainties in the evolution of the SBHBs with the smallest (sub- 
parsec) orbital separations; such binaries are also the most inter- 
esting because of their proximity to gravitational wave regime 
and coalescence. Further improvements in the treatment of gas 
thermodynamics can be made thus reducing one of the two 
dominant sources of uncertainty (the other one being the initial 
conditions). This can be achieved by adopting a more physi- 
cally motivated form for gas cooling and also gas heating by 
introducing the sources of ionization associated with the accret- 
ing SBHs. 

Given the challenges associated with the detection of SB- 
HBs and recoiling SBHs, predictions of observational signa- 
tures based on existing hydrodynamical simulations would be 
very valuable. This is currently not attainable because hydrody- 
namical models do not include radiative transfer prescriptions 
that allow such predictions to be made (the exception are the ac- 
cretion rate curves, which are also subject to the uncertainties 
described in the previous paragraph). Indeed, a complete treat- 
ment of radiative transfer within a hydrodynamical simulation 
of SBHs and non-zero metallicity gas is currently not feasible. 
Nevertheless, a significant step toward this goal can be made 
with the help of photoionization codes that include full radiative 
transfer treatment of gas composed of multiple atomic species 
and can make a variety of predictions in the form of character- 
istic emission lines and continua, associated optical depths, and 
a plethora of other information. Photoionization codes can be 
used to characterize radiative transfer within every individual 
parcel of gas in a simulation, given the physical properties at the 
location of a parcel (number of ionizing photons, temperature, 
and density of the gas). Because the range of physical proper- 
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Figure 5: Ha emissivity of the gas as a function of radius, at 9.4 years for a 
co-rotating model witli solar metallicity gas. The emissivity of each gas cell 
plotted in the figure is weighted by the density of the gas at that position, so 
that comparison with the parametric emissivity model (plotted as a solid line 
with arbitrary normalization) can be mad e. The color legend is the same as in 
previous two figures. Figure adapted from lBogdanovic et alj j2008l) . 



ties in a simulation can be predicted with some certainty (based 
on the properties of nuclear regions in galaxies hosting AGNs or 
from test simulations), the response of the gas to photoioniza- 
tion precomputed and stored in a database. Such a database can 
be used for re ference during sim ulations (for an example of this 
approach see lBogdanovic et al llioos), or in some cases in the 
post-processing of simulation data (IBogdanovic et al.Ll2004l) . in 
order to calculate observational signatures. The main effort as- 
sociated with this approach is the calculation of well behaved 
heating, cooling, and emissivity maps (within the limits of the 
specific photoionization code). The computational cost associ- 
ated with the addition of the photoionization data grid to the hy- 
drodynamical code is modest, as the majority of the additional 
operations are interpolations between the grid points. While 
any such method must acknowledge that detailed modeling of 
emission line and other spectral properties is currently still out 
of reach, it still offers the best means of examining the most 
characteristic signatures of sub-parsec SBHBs and recoiling 
SBHs. 

In the following section we summarize the results from such 
a study of observational signatures of SBHBs interacting with 
gas. 



3.1. Investigations of observational signatures using multi- 
species radiative transfer 

We implemented an approximate multi-species radiative 
transfer scheme in a modified version o f the smoothed parti- 



cle hydrodynamics (SPH) code Gadget 1 (ISpringel et al.Ll2001h 
in order to study the observational sig natures of sub-parsec bi- 
naries interacting with gas (Bogdanov ic et al The gas 
is assumed to have solar metallicity and its physical proper- 
ties were characterized by calculating heating, cooling, and ra- 
diative processes as an integral part of hydrod ynamical sim- 
ulations. Using photoionization code Cloudy (iFerland et al. 



Figure 6: Effective accretion rate on the primary {solid, red line) and secondary 
{dashed, blue line) black hole calculated from co-rotating model. Note that 
in this figure only we show data from a lower resolution test simulation (20k 
particles) which ran over a longer period of time (~ 50 yr) than the base- 
line look simulation (~ 35 yr) but is in all other aspects equivalent to it. 
The accretion rates can be translated into UV/X-ray luminosity by assuming 
1 Mq yr"' » 10^^ erg s"' in the UV/X-ray band. The arrows mark the times of 
the pe ricentric passages of the binary. Figure adapted from , Bogdanovic et al] 
<2008l). 



1998h . we constructed a grid of photoionization maps for a 
range of parameters, determined from preliminary simulations 
without cooling. The maps are calculated in the parameter 
space of density and temperature of the gas and intensity of the 
ionizing radiation. The range of parameter values for which 
the maps were computed is as follows: 10^ cm"-' < n < 
lO'^cm ^ 2000 K < T < 10*^ K, and erg em 's ' < 7 < 
lO'^ erg cm"^ s"'. Gas at higher temperatures than 10^ K is 
also found in our simulations; in such cases we calculate heat- 
ing and cooling rates and other properties of interest by linearly 
extrapolating the grid values. In the SPH simulations we set 
the lower threshold for the gas temperature to 100 K and as- 
sign no upper threshold. Based on these results we calculated 
the accretion-powered continuum and Hot light curves, as well 
as the Ho- emission line profiles emerging from the inner parts 
of a gas disk in the presence of a binary, on a scale of < 0.1 
pc. We selected the Balmer series Ho- line (/irest = 6563 A) 
because it reflects the kinemati cs in the broad line region of 
AGN (e.g.. lSulentic et akluOO O). Hence, by mapping the kine- 



matics in the region of sub-parsec size, broad He line can in 
principle provide important information about SBHBs interact- 
ing with gafl Moreover, Hff is the most prominent broad line 
in the optical spectrum, least contaminated by neighboring nar- 
row emission lines. The Ha line in accretion-powered sources 
is thought to be powered by illumination of a disk by an AGN 
and the radi ative transfer associated w i th it i s relatively well 
understood dComn-Souffrin & Dumoiil Il989l) . That said, the 
method can be extended to include other transitions of inter- 



^We do not model the narrow component of Hff line which is thought to 
originate at larger distances. 
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Figure 7: Trailed speclrogram (left) and black hole projected velocity curves 
(righl) plotted for the co-rotating model. Trailed spectrogram is a logarithmic 
gray scale map of the Ha intensity against wavelength and observed velocity. 
Darker shades mark higher intensity. Notice the low relative intensity of the 
profiles before the start of accretion. The velocity curve panel shows the orbital 
velocities of the primary (dashed, red line) and secondary (solid, blue line) 
black holes projected along the line of sight to the observer The velocity curves 
are skewed because of the non-zero eccentricity of the orbit. Characteristic 
features associated with the peiicentric passages of the binary in years 8 an d 23 
are marked on both figures. Figure adapted from lBogdanovic et alj )2008l) . 



est, given that calculation of the photoionization data grid with 
Cloudy includes a wide range of atomic species from hydrogen 
to zinc. 

3.1.1. Physical properties of the gas: 

In the remainder of this section we focus on selected results 
from the simulation of the binary co-rotating with the gas disk 
and refer the reader to the original publication for discussion of 
the counter-rotating case. The simulation follows the evolution 
of a SBHB over 35 years in total and comprises 10^ SPH gas 
particles (100k hereafter). The black hole mass ratio is 10:1, 
with the primary having a nominal mass of Mi = 10** M©. The 
secondary SBH is on an eccentric orbit around the primary with 
e = 0.7, and the initial orbital period of the binary is close to 16 
years. The gas disk is initially associated with the primary black 
hole only and its mass and outer radius are initially lO'* and 
2 X lO^r^ {rg = GMi/c^). The median value of the initial disk 
density is ~ lO'^ cm~^. The geometry of the system and the 
temperature of the gas are illustrated in Figure 3 at f = 9.4 yr 
after the beginning of the simulation. The primary SBH is ini- 
tially located at {x, y, z) = (0, 0, 0) and exhibits modest orbital 
motion during the simulation. 

The gas in the simulation is heated by shocks that arise from 
interactions with the secondary SBH and by radiation once the 
accretion onto the black holes begins. As a result, a portion 
of the gas is blown out of the plane of the disk and forms a 
hot halo that envelopes the disk. The gas quickly departs from 
the uniform initial temperature of 2000 K and assumes a wide 
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Figure 8: Sequence of Ho' emission-line profiles selected from a model where 
the disk is co-rotating with respect to the orbit of the SBHB. The intrinsic inten- 
sity of profiles is plotted against wavelength. The first 4 profiles in the sequence 
are multiplied by a factor of 10. so that they can be represented on the same in- 
tensity scale with the other profiles. The con'esponding time from the beginning 
of the simulation is given next to each profile. The observer is located at infin- 
ity, in the positive xz-plnne at an a ngle of i = 30° wi t h resp ect to the axis of the 
binary orbit. Figure adapted from lBogdanovic et all i2008h . 



range of temperatures, as shown in Figure 4. Figure 4 further 
illustrates that multiple temperature components can be present 
at a single radius and shows the temperature distribution of the 
SPH particles. The hot component of the gas spends a signifi- 
cant amount of time in the temperature range 10"* - 10** K and in 
this regime radiative cooling is dominated by free-free emission 
and recombination. The highest value of the gas temperature in 
the simulation, T~ 10'^ K, is reached after the shock is formed 
by the secondary. On a time scale of months, the temperature 
of the gas falls below 10'° K due to the combined effects of 
radiative cooling and adiabatic expansion. On the other hand, 
the cold component of the gas, confined to the higher density 
spiral arm has a temperature in the range 10^ - 10^ K, close to 
the initial value, and cools largely through recombination radi- 
ation. At the end of the simulation the gas density in the disk 
is in the range 10^ - 10'"* cm"^. The median density in the ten- 
uous, photoionized halo is 10^ - 10** cm and the minimum 
density can be as low as 10'' cm"^. 

The interaction of the binary with the gas also affects the 
emissivity of the gas and consequently its emission signatures. 
The most important consideration in the calculation of the emit- 
ted Ha light is the efficiency with which the gas reprocesses the 
incident ionizing radiation and re-emits it in the form of Ho; 
photons. This efficiency is commonly characterized by the sur- 
face emissivity of the gas. The emissivity of the photoionized 
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gas depends on numerous physical parameters. Locally, it de- 
pends on the physical properties of the gas. Globally, its spatial 
distribution depends on the structure of the accretion disk. We 
assess the Ha emissivity of every gas particle in the simulation 
using photoionization data grids precomputed with Cloudy. In 
Figure 5 we compare the emissivity determined from our pho- 
toionization calculations with the parametric model used for 
AGNs in which the emissiv ity is described as a function of ra - 
dius, e ~ ^ * with q ^ 3 jComn-Souffrin & DumonI Il989h . 
We find a qualitative agreement between the two albeit, with a 
significant amount of scatter This is not surprising given the 
variations in surface density of the perturbed gas disk in our 
binary model. 

3.1.2. Characteristic observational signatures: 

We further find that X-ray outbursts (powered by accretion 
onto the two AGN and shocks) should occur during pericen- 
tric passages of a coplanar binary as long as the nuclear re- 
gion is not devoid of gas. During pericentric passages the bi- 
nary accretion rate reaches ~ 30 yr"' (Figure 6), compa- 
rable to the Eddington rate for this system, while during the 
remainder of the orbital period the total accretion rate aver- 
ages at the level of few Mq yr"'. This implies luminosity of 
Le ~ 1.5 X lO^'^Ms erg s"' during outbursts and an order of 
magnitude lower luminosity on average. Despite the 10:1 mass 
ratio the two SBHs exhibit comparable accretion rates and at 
some times the accretion rate onto the secondary exceeds that 
of the primary SBH. This inversion happens the lower mass ob- 
ject has a smaller relative velocity with respect to the gas than 
the primary and c an result in a higher relative luminosity of a 
lower mass object ("Artvmow icz & Lubow , 1996 : Gould & Rix , 
I2OOO; Hayasaki et al., 2007). During this same period the inte- 
grated (bolometric) luminosity of the disk is at a mean level of 
~ lO'^'^ erg s and thus, comparable to that of the photoioniza- 
tion sources. 

At the mean level of accretion inferred from the simulation, 
the gaseous disk should be depleted after only lO^* - 10"* yr In 
order for the outburst activity to last over longer time scales, 
the reservoir of gas in the nuclear region needs to be contin- 
ually replenished. It is also plausible, however, that repeated 
collisions of the secondary black hole with the disk will com- 
pletely disrupt it and turn it into a spherical halo of hot gas. In 
such a case, the accretion rate would be relatively smooth and 
uniform, and no outburst would be evident. A calculation fol- 
lowing the X-ray light curve variability over a large number of 
orbits is necessary in order to confirm that the periodicity is a 
long-lived signature of the binary. Although currently not pos- 
sible with the SPH method used here, simulations of the long 
term evolution of the binary together with hydrodynamics and 
radiative transfer may be achieved in the future. 

In addition to the recurrent outbursts in the X-ray light curve 
the signature of a binary is potentially discernible in the Ha 
light, and more specifically. Ha emission line profiles. After the 
start of accretion, the Ha luminosity reaches 10-^''- 10"*" ergs"', 
observable out to the distance of the Virgo Cluster (16 Mpc) 
and possibly up to the distance of the Coma cluster (100 Mpc). 
The changes in the profile sequence are easiest to discern in the 



trailed spectrogram (shown in the left panel of Figure 7) which 
represents a 2D map of the Ha intensity against observed veloc- 
ity (or wavelength). The Doppler boosting of the blue side of 
the line and gravitational redshift of the red wing are noticeable 
during the pericentric passages of the binary and, in principle, 
can serve as indicators of the orbital period. The width of the 
Ho- profile increases after the pericentric passages of the sys- 
tem, reflecting the inflow of gas towards the primary. Also, the 
widening of the profile is asymmetric and slightly shifted to- 
wards the red with respect to the pre-pericentric sequence of 
profiles. This shift is a signature of the motion of the accretion 
disk which follows the primary that is receding from the ob- 
server In Figure 7 we compare the trailed spectrogram with the 
velocity curves of the two binary components projected onto 
the line of sight to the observer and highlight the instances of 
pericentric passages. The variations in the profile intensities 
correspond to the features in the velocity curve of the secondary 
in years 8 and 23, when it is moving towards the observer with 
the highest projected velocity (i.e., the binary is in quadrature). 
Hence, we find the signatures of the velocity curves of both 
black holes in the emission line profiles: if such features could 
be discerned in the observed Ha emission-line profiles, they 
would signal the presence of a binary and potentially allow for 
determination of its mass ratio. 

The sequence of profiles shown in Figure 8 can be obtained 
if the trailed spectrogram is sampled at selected times. Initially, 
from the unperturbed disk we observe double-peaked emission- 
line profiles but they gradually depart from this shape as the 
perturbation propagates through the gaseous disk. Our calcu- 
lated profiles account only for the broad component of the Ha 
line; the narrow component is not taken into account. The nar- 
row component of Ha is present in the observed (real) profiles 
and is thought to originate in gas outside of the nuclear region. 
Its presence adds a level of complexity to the analysis and mod- 
eling of the observed Ha emission-line profiles, as the narrow 
component of the line needs to be removed before the kinemat- 
ics and geometry of the gas in the nuclear region can be assessed 
from the broad component of the profile 

The irregular and variable profiles in Figure 8 are not unlike 
those seen in a small percentage of AGN and in SBHB can- 
didates found in SDSS. In practice however, going from the 
detection of a binary candidate to a positive identification of a 
SBHB based on only a few epochs of observations of irregular 
broad Balmer lines is quite a challenging task. Our study in- 
dicates that the variability of broad emission-line profiles may 
be subtle and reflected in the extended wings of the profiles. In 
that sense an expectation that the orbital motion of the binary 
can be inferred from the velocity shift of the profile peaks may 
be simplistic and it may yield unreliable results. 



^^For a more detailed discussion of phenomena that may affect the shape of 
the Ha emission l ine profiles as wel l as th e effect of approximations made in 
the calculation see lBogdanovic et alj )2008l) . 



7 



4. Conclusions 

Understanding the population of SBHBs near coalescence 
in the local universe can provide insight into the population 
at high redshif t where theoretical m odels suggest most merg- 
ers take place (iSesana et all l2007bl) . In practice however, the 
search for close binaries and their post-coalescence counter- 
parts based on their EM signatures has proven to be a difficult 
task. There is only a handful of observed bound SBHB candi- 
dates and no confirmed cases at all. Thus, our understanding of 
the accretion physics and associated observational signatures is 
poor. Given the challenges associated with detection of SBHBs 
and recoiling SBHs it would be valuable to make predictions 
of related observational signatures using hydrodynamical sim- 
ulations. While some steps have been taken in that direction, 
more development is needed in the future in order to calculate 
multi-band observational signatures of these objects. 

Our pilot study implies that the period and perhaps the mass 
ratio of the binary can be measured from well-sampled, long- 
term X-ray and optical light curves and Ha profile sequences 
that have been followed for at least a few revolutions of the bi- 
nary. In order to achieve an efficient SBHB search and avoid 
time-consuming observational foUowups of non-SBHB candi- 
dates (false positives) it is important to utilize all available com- 
plementary signatures that may help to break this degeneracy. 
Given the location of selected SBHB and recoiling SBH candi- 
dates on the sky, the best prospect for detection is with radio in- 
terferometers, if both black holes have associated radio sources. 
More immediate tests of existence and coalescence of binaries 
may be possible in the future in synergy with GW observations, 
once the Laser Interferometer Space Antenna become opera- 
tional. 

We would like to thank the anonymous referee for insightful 
comments and suggestions. TB thanks the UMCP- Astronomy 
Center for Theory and Computation Prize Fellowship program 
for support. 
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